Quantum search by partial adiabatic evolution 

Ying-Yu Zhangj and Song-Feng Lu 
School of Computer Science, Huazhong University of Science and Technology, Wuhan, China. 

(Dated: July 12, 2010) 

A quantum search algorithm based on the partial adiabatic evolution,!] is provided. We calculate 
its time complexity by studying the Hamiltonian in a two-dimensional Hilbert space. It is found 
that the algorithm improves the time complexity, which is 0{^J N/M), of the local adiabatic search 
algorithmg], to 0{Vn/M). 



PACS numbers: 03.67.Lx 
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I. INTRODUCTION 

Quantum adiabatic computation has attracted a lot of 
attention in the past decades, such as [ll. l3l-[ig| . since it 
was proposed by Farhi et al. [ll| . In (5| , an adiabatic algo- 
rithm was proposed to solve the Deutsch-Jozsa problem. 
The algorithm took an exponential time which provided 
only a quadratic speed up over the best classical algo- 
rithm for the problem. A modified algorithm proposed 
by Wei and Ying |12j | improved the performance to con- 
stant time. In [6|, quantum adiabatic computation was 
proved to be polynomially equivalent to the quantum cir- 
cuit model. The proof showed that adiabatic quantum 
computation using Hamiltonians with long-range five- 
or three-body interactions, or nearest-neighbor two-body 
interactions with six-state particles, could efficiently sim- 
ulate the circuit model. This results was soon modified to 
qubits with two-body interactions [ij]. A simpler proof 
of the Equivalence was presented in [8| . In [1J| , quantum 
adiabatic computation was applied to solve random in- 
stances of NP-complete problems. A research outline of 
its application to solve NP-complete problems can also 
be found in the paper. 

A typical quantum adiabatic algorithm starts with the 
ground state of the initial Hamiltonian Hi, and evolves 
slowly to the ground state of the final Hamiltonian Hf . 
The system that implements the algorithm uses the time- 
dependent Hamiltonian 



H{s(t))^{l-s{t))H.,+s{t)H}. 



(1) 



search problem [l5|, and designed a quantum search al- 
gorithm based on local adiabatic evolution. In the algo- 
rithm, Hi and Hf are specified as 



and 



where 



H, = l-\^){^, 



Hf = l^\P){P\, 



N-l 



N 



(4) 



(5) 



(6) 



The running time (evolution time) is essentially deter- 
mined by [l| 



and 1/3) is an equal superposition of all marked states. 
The algorithm can find one marked item in a running 
time of order ^J N /M . 

In [1] , Tulsi proposed a partial adiabatic evolution with 
Hf a one-dimensional projector Hamiltonian. It was also 
checked in the paper that Roland and Cerf 's results (of 
the case that there is only one marked state) can be ob- 
tained as a special case of the partial adiabatic evolu- 
tion. In this paper, we give a specified quantum search 
algorithm based on the partial adiabatic evolution, and 
show that the algorithm provides a better time complex- 
ity, which is 0{Vn/M), than that, which is 0{^N/M), 
of the local adiabatic search algorithm. 

This paper is organized as follows. In section 2, we 
specify the partial adiabatic search algorithm. In section 
3, we calculate its time complexity. We conclude the 
paper in section 4. 



T>e(-^) 
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Jmin 



where 



(2) 



g™„ = min[E{l, s) - E{Q, s)], (0 < s < 1). (3) 

E{0,s) and E{l,s) are the two lowest eigen values of 
H{s). Roland and Cerf [2] considered the unstructured 
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II. PARTIAL ADIABATIC SEARCH 
ALGORITHM 

In this section, we specify the partial adiabatic evolu- 
tion [l| as a search algorithm. For convenience, we use 
different notations. The search algorithm executes the 
four steps below: 

(1) The initial state is prepared to be |^). 

(2) At t=0, the Hamiltonian is suddenly changed to 
H{s^) without disturbing the state j'l'). 

(3) The Hamiltonian evolves from H{s^) to H{s^) lin- 
early in time over duration T' . 



(4) Measure the state of the system. 
Repeat these four steps until we find a marked state. 
For the algorithm, s~ and s+ are specified as s^ = ^ — 
— W= and s+ = i + —ttt- Because the algorithm evolves 

2\/J\ ^ 2\/J\ 

adiabatically only within a small time interval [s~ , s"*"], it 
is called a partial adiabatic search algorithm. After step 
2, the system that implements the algorithm will be still 
in the state |\1/). That is, the system state will be the 
ground state of H{s~) with probability |(*|£;(0, s~))p. 
The adiabatic theorem \W\ guarantees that it will be the 
ground state of H{s+) with probability |(^|S(0, s"))p 
after step 3. Measuring the state of the system will give 
the ground state of Hf = 1 — |/3)(/3| with probability 
P = |(*|£;(0,s-))|2 X |(/3|i;(0,s+))|2. We call P one- 
round success probability, and accordingly T' one-round 
evolution time. The overall time complexity(evolution 
time) of the algorithm is T = T'/P. Here, T' = uj/g: 
where cj = s+ — s^ = -4=. 

VN 



III. TIME COMPLEXITY 
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Let 



minimum energy gap 



|"^ = ^^Ei-) 



VN - M 



xg:S 






(7) 
(8) 



x£S 



where S is the set of the marked states, and M is the 
number of the marked states. Throughout of this paper, 
we suppose S is not empty. This means that \ < M < N . 
The state |^) can be rewritten as [17 1 



I*) 



N - M. 



N 



a 



(9) 



It is easily to check that any state orthogonal to |a) and 
1/3) is a eigen state of H{s) and the corresponding eigen 
value is 1 which is (N-2) times degenerated. This means 
that E{k, s){k — 0, 1) and are in the subspace spanned by 
\a) and |/3). Besides, all the following calculations involve 
states only in the two-dimensional subspaces panned by 
I a) and |/3). So, we can work in the subspace, instead of 
working in the N-dimensional Hilbert space. 
Let the eigen spectrum of H(s) be 



H{s)\E{k,s))=E{k,s)\E{k,s)), 



(10) 



where E{k, s) and \E{k, s)) are the k-level eigen value 
and eigen state of H{s) respectively. Throughout this 
paper, we only consider the two lowest eigenstates and 
eigenvalues, i.e. k=0,l. Left multiplying {a\ to Ea. ([TU| . 
we get 

{a\H{s)\E{k, s)) = E{k, s){a\E{k, s)). (11) 



Substituting the Eq.([T]) into the left side of Eq. lfTTj) . we 
also get 

{a\H{s)\E{k,s)) = 
{a\E{k, s)) - (1 - .s){a\^)mE{k, s)), (12) 

with (a|/3) = 0. Combining Eas. (fTT|) with (IT2]) . we obtain 

^^l^^'^' '^^ ^ l~Eik,s) ' (^^) 

when 1 — E{k, s) ^ 0. Similarly, 

(WM))^ "-"<^l''>gf'^-» . (14) 

1 — s — E[k, s) 

when 1 - s - E{k, s) ^ 0. Thus, with Eqs.dlSl) and (flil) . 
we get 

{^\E{k,s)) = (vl/|a)(a|i?(fc,s)) + (vl/|/3)(/3|i?(fc, s)) 
(1 - s)A (1 - s)B 



= ( 



\-E{k,s) l-s-E{k,s) 



){^\Eik,s)){l5) 



where A == |(*|a)|2 = ^^ and B = |(*|/3)|2 ^ §. Re- 



arranging Eq. ljlSp . we finally obtain the secular function 
for H(s) 



E^{k,s)-E{k,s) + s{l-s)A^O. (16) 



So, 



and 



EiO,s),Ei 



M) . i^v^lHifliEill, (1 



7) 



= E{1, s) - E{0, s) = Vl - 4s(l - s)A. (18) 



The minimum energy gap gmin = \f\ — A — ^J M/N is 
obtained for s = 1/2. The one-round evolution time is 

T = y/N/M. 

B. one- round success probability 

In this subsection, we calculate the one-round suc- 
cess probability. Substituting Egs. p^ and P^ into 

\{a\E{k,s))\^ + \{(3\E{k,s))\^ = l(fc = 0,1), we have 

, (1-5)2^ (l-s)'^B ,,, , , ,,,, 

^ n-EiLsw + (i-L4,.))2 )'<^'^(^'-)^' - ' 

(19) 
whenl-£;(yfc,s) y^ and l-s--B(fc, s) ^ 0. If I-st^O, 
this immediately gives 

mE{k,m' 



(-'- ^' ((l-B(fc,s))^ + (l-s-B(fe,s))^) 



Substituting Eq.dlOl) into Eq.dTH), we get 



B 



(20) 



i 1-E{k.s) ) A + B 



Because 



A 



{l-s-f = \{l + ^f<l, (22) 



<4, (23) 



4(iV - M)/N 



{l-E{Q,s-)Y n + /i+M_ jvTp 



and 



B 



AM 



(1-s -EiO,s ))2 (1+ /l + M„M) 



< 



4Af 



(1 + VA'/) 



/n2 



we obtain 



using Eq. pO)) . 
Because 

l-s^ 



|(*|£;(o,.-))p> 



^(0,g+) _ 7F + \/ 



M+l _ M 



1 - E{0, 



M+l _ M 
N "W 



< 



-1 



N 



(24) 



<4, (25) 



(26) 



^^^-4<a/^ = V2^, (27) 



we obtain 

|(/?|i?(0,.+ ))p> 



B 



^ >i (28) 



2BA + B 2A + 1 3 



using Eq. (PT|) . 

This provides a lower bound of the one-round success 
probabihty 

P=mE{0,l ^))|2x|(/3|£;(0 i + ^))|2>-^. 

2 2^/N 2 2V^ 24 

(29) 

So, the overall time complexity of the partial adiabatic 

search algorithm is T = T'/P = 0{Vn/M). 



IV. CONCLUSION 

We have provided a quantum search algorithm based 
on the partial adiabatic evolution. As we have seen that 
the minimum energy gap along with the one-round suc- 
cess probability determine the overall time complexity 
of the algorithm. They are calculated by studying the 
time-dependent Hamiltonian H{s) in a two-dimensional 
Hilbert space. The overall time complexity is 0{\/N/A4), 
which provides a speed up of y/M over the local adiabatic 
search algorithm. This is resulted from the fact that the 
one-round success probability is bounded from below by 
a constant. 
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